Signature of the FFLO phase in the collective modes of a trapped ultracold Fermi gas 
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We study theoretically the collective modes of a two-component Fermi gas with attractive inter- 
actions in a quasi-one-dimensional harmonic trap. We focus on an imbalanced gas in the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phase. Using a mean-field theory, we study the response of the 
ground state to time-dependent potentials. For potentials with short wavelengths, we find dramatic 
signatures in the large-scale response of the gas which are characteristic of the FFLO phase. This 
response provides an effective way to detect the FFLO state in experiments. 



With the advent of ultracold atomic gases new states 
of matter are becoming experimentally accessible. Im- 
portant developments have included the achievement 
of superfluidity in two-component Fermi gases [H, and 
the study of the effects of density imbalance in these 
systems [2|. These developments open up the prospect of 
the study of unconventional superfluid phases, in partic- 
ular the inhomogeneous superfluid state of Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) [1]. This state has been 
under theoretical investigation since the 1960s, playing 
a central role in our understanding of superfluidity of 
fermions both within condensed matter physics and in 
particle physics yet to date no unambiguous demon- 
stration of the FFLO state has been achieved. 

While theory predicts the appearance of the FFLO 
phase in a bulk (3D) gas of ultracold fermions j3] the 
region of parameter space that it occupies is expected 
to be very small. It has been shown that this phase is 
greatly stabilised in (quasi-) ID geometries [1, 0, B 0- 
Such geometries can readily be achieved in atomic gases 
by optical confinement [lo| . so experiments on ultracold 
Fermi gases will soon be in the regime where the FFLO 
phase can be studied. In view of this, it is important and 
timely to ask: what are the observable properties of the 
FFLO phase? 

Existing proposals for the detection of the FFLO phase 
include the probing of noise correlations in an expan- 
sion image following release of the trapQ, and the use 
of RF spectroscopy to excite atoms out of the gas into 
other states ll| . In this Letter we show that characteris- 
tic signatures of the spatial inhomogeneity of the FFLO 
phase can be found in the collective mode response of the 
trapped gas. The method directly probes the intrinsic 
inhomogeneity of the FFLO phase - by the use of a per- 
turbation by an potential with a short-wavelength spatial 
periodicity. An important feature of this method is that 
the signatures of the microscopic nature of the phase ap- 
pear in the response on very large length scales (the scale 
of the atomic cloud) . Consequently, the measurements of 
the response do not require high spatial resolution, and 
could be performed in situ without expansion of the gas. 

We study a model of two-component fermions in ID 
with attractive contact interactions, and unequal densi- 
ties. For a homogeneous system, the exact groundstates 



are known from the Bethe ansatz 0, ■ Depending 
on the densities of the two components there are three 
phases that appear at T = 0: the unpolarised superfluid 
state; a fully polarised (therefore non-interacting) Fermi 
gas; and a partially polarised phase 0], associated with 
the FFLO state. This is a state that, within mean-field 
theory, has a local superconducting gap that oscillates in 
space. As has been discussed in detail by Liu et al.[l^, 
mean-field theory provides an accurate description of the 
exact phase diagram for the ID system, at least within 
the weak-coupling BCS regime. We have extended the 
mean-field theory to investigate the linear response and 
collective modes of both the trapped and untrapped im- 
balanced Fermi gas. As we shall discuss in detail later, 
this approximate theory correctly captures all qualita- 
tive features of the collective modes that are important 
for our purposes. 

To effect numerical calculations, we study a dis- 
cretized version of the problem (the attractive Hubbard 
model), for which the mean-field Bogoliubov-de Gennes 

(BdG) Hamiltonian is H = —JJ2i a ('^l+i cr'^i,'^ + 



where c,-t! are 



fermionic operators for species a =t, J. on site i, Wi a 



jia (with y^^* the external potential 



and /Ict the chemical potentials), and = U (ci^iCi^i) 
is the local superfluid gap. J is the hopping parameter 
and U the on-site interaction strength (?7 < is assumed 
throughout). We derive the linear response to external 
time-dependent perturbations, SWi^ait), by supplement- 
ing the (self-consistent) solutions of the above BdG equa- 
tions with the random phase approximation (RPA) [ij]. 
Divergences in the response appear at the frequencies 
of the collective modes of the system. The results we 
present are at sufficiently low particle density to be rep- 
resentative of the continuum limit. In the mapping to 
the continuum, we relate site i to position x via x ~ ia 
and the mass is m = TTa^' interaction strength 7, 
defined as the ratio of the interaction energy density to 
the kinetic energy density, is given by 7 = — ™f p'' [isj 
where gid ~ — , p is the total density of particles and 

We discuss first the response of a homogeneous sys- 
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FIG. 1: Response [I^ of a homogeneous system in the FFLO 
phase to a spin-asymmetric periodic potential ([T]). The area 
of the circle is proportional to the amplitude of the response. 
The polarisation p = ^^^j^^^j = 0.15, and the interactions 
strength is 7 = 1.5. Two gapless sound modes are seen to 
emerge around fc ~ and ~ fc* = 2(fcF^ — kpi)- The 
simulation was done on 270 lattice sites. 



tern, with = and periodic boundary conditions. 

We focus on the FFLO phase, with unequal average par- 
ticle densities, for which the self-consistent mean-field 
solution has an oscillating gap A (a;) with wavelength 
Aa = 27r/(fci?.| — fci^^^i), where kp^^r are the Fermi wavevec- 
tors for the non-interacting imbalanced gas. We study 
the response to periodic perturbing potentials 



6Wa {x,t) = Va sin kx cos ujt 



(1) 



We refer to the case ~ Vi ~ Vq as "spin-symmetric" 
and V| = — = Vo SiS "spin-asymmetric" excitation. 

The response of an imbalanced gas in the FFLO phase 
to a spin-asymmetric modulation is shown in Fig. [1] 
At low frequency and long wavelength there appear two 
distinct sound modes with different velocities [3, 17 1. 
Within mean-field theory, the appearance of these two 
gapless modes arises from the fact that the FFLO phase 
breaks both gauge symmetry and translational symme- 
try. (We later discuss these properties beyond mean field 
theory.) An analysis of the two low-frequency modes 
shows that they are of mixed character, involving spa- 
tial oscillations of both density and spin-density. This 
is expected from the breaking of time reversal symmetry 
by the imbalance 17| , and is consistent with the results 
(in Fig. [1]) that both modes can be excited by a purely 
spin-asymmetric perturbation, V| = — 

In Fig. [1] it is clear that linear sound modes emerge also 
around the point k = k* = 2{kp-^ — kp^). Within mean- 
field theory, this too can be understood as a consequence 
of the broken translational symmetry of the FFLO phase. 
This leads to a Brillouin zone structure for the collective 
modes, characterised by a reciprocal lattice vector of size 
k* . Note that this is twice the value that one would ex- 
pect from the translational periodicity of the gap, Aa- To 
understand this, observe that, while the gap A has pe- 
riod Aa, the densities Pi{x) and Pi{x) have periods Aa/2, 
and the mean-field ground state has an exact symmetry 
under the transformation x x + Aa/2, A —> — A. Ap- 
plying this symmetry to the RPA response calculation 



leads to a generalised Bloch theorem for the collective 
modes. The wavevector k is found to be conserved mod- 
ulo (27r)/(AA/2) = k*, consistent with the response in 
Fig. [TJ For a spin-symmetric perturbation a similar qual- 
itative response is found. However, the response close to 
fc = fc* is smaller than for the spin-asymmetric pertur- 
bation. We account for this by the fact that the periodic 
density of excess majority particles manifests itself more 
strongly in the difference of the two densities than in the 
sum of the two densities. We now turn to discuss the 
manifestations of this response for a trapped gas. 

We have studied the attractive Fermi gas in a har- 
monic trap, with V^^^{x) = ^mwgX^ and with unequal 
particles numbers iV| 7^ 7V|. According to Orso [6| a 
ID imbalanced trapped Fermi gas can only be in one 
of two configurations: (a) A partially polarised phase in 
the middle, with a fully paired phase towards the edge; 
(b) a partially polarised phase in the middle, with a 
fully polarised phase towards the edge. These configu- 
rations appear within the self-consistent BdG mean-field 
theory [isj. An example of the case (a) is given in Fig. [21 
and of the case (b) in Fig. [31 where lengths are measured 
in units of N^^^Uho (with N = + Ni the total number 
of fermions, and aho = \/hJrnuj^ the oscillator length of 
the trap). The configurations in Figs. [2l and [3l are in the 
weak coupling ECS regime, having values of 7 for which 
BdG gives qualitatively the right density profiles 
We have studied the response of the trapped imbalanced 
Fermi gas in both these regimes. 

Perturbations with wavelength much longer than the 
size of the system couple via potentials Vaix) that are 
proportional to x. We find that the response to such 
potentials is well described by two sharp modes, involv- 
ing the dipolar oscillations of the density and spin. One 
of these modes - the "Kohn mode" - involves the in- 
phase motion of both spin components, and has fre- 
quency Lo ~ independent of the state of the system. 
This exact result holds also within RPAflsj. and is re- 
covered to high accuracy in our calculations, showing 
that discretisation effects are minimal. The other mode - 
the "spin-dipole" mode ~- involves the relative motion of 
the two spin species. For attractive interactions the fre- 
quency lies slightly above wo (at lo^a = I.SSwo and 
for Figs. [21 and [21 respectively). While the frequency of 
the spin-dipole mode does depend on the state of the 
system, it is does not show strong features of the pres- 
ence of the FFLO phase. Similarly, the breathing modes 
are insensitive to the microscopic nature of the atomic 
gas0. 

The most interesting features arise when one excites 
the gas by time dependent potentials H]) with short wave- 
lengths. We find characteristic signatures of the micro- 
scopic nature of the trapped gas in the response of spin- 
dipole mode " i.e. on length scales much larger than the 
wavelength of the applied perturbation. 27| In Figs [H^b) 
and [2Kb) we show the amplitude of the response of the 
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FIG. 2: (Colour online.) Configuration with fully paired 

phase towards the edge of the trap, (a): densities and the 

value of the superfluid gap A. (b): response of the spin- 

dipole mode to excitations of different wavelengths (arbitrary 

units). Here p — 0.048, 7 — 0.93 (measured in centre), num- 

_o 1 

ber of particles A'^ = 290, lattice spacing a = 3.3-10 N^aho- 
The perturbing potential has a fixed amplitude, while the 
wavelength is varied. 
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FIG. 3; (Colour online.) Configuration with fully polarised 
phase towards the edge. Herep = 0.25, 7 = 1.5, a/{N^aho) ~ 
4.7 ■ 10"^, iV = 143 . Otherwise the same as fig.H 



spin-dipole mode to spin-symmetric and spin-asymmetric 
periodic perturbations ([T]) as a function of wavelength 
A = 27r/fc. 

For wavelengths A that are short compared to the 
size of the cloud, A^^/^Oho, but large compared to 
Aa, we find an oscillatory response as a function 
of A in both Figs. EJb) and E^b). To understand 
this behaviour, we note that the spin-dipole mode 
can be excited when the two species experience a 
net difference in their acceleration. This difference 
is a{k) / [p| (a;)/iV-|- =F Pi{x)/Ni]cos{kx)dx where 

the =F signs denote spin-symmetric/spin-asymmetric ex- 
citations. The difference in acceleration oscillates with fc, 
in a manner that is largely controlled by the sizes of the 
clouds. For case (b) we assume, for simplicity, parabolic 
density distributions of size Xa for the two spin compo- 
nents a [2I]. Defining r = ^^^^^ and S = , the 
condition for a maximum is sin kr cos kS = for spin- 
asymmetric excitations, and cos kr sin kS = for spin- 
symmetric excitations. For S r it follows that for 
spin-asymmetric (spin-symmetric) potentials maxima in 
the response come from the sin kr (cos kr) term and the 
antinodes in the envelope come from the cos kS (sin kS) 
term. Case (a) can be described by assuming a parabolic 
density distribution for the average density but where 



there is a constant difference between the two densities 
in the partially polarised region of size rp. For spin- 
symmetric and spin-asymmetric excitations the condition 
for a peak is given by cos krp = 0. The modulation of the 
peaks by this envelope structure allows us to distinguish 
the two imbalanced trap configurations (a) for which the 
two clouds have equal size x-\ = xi and there is no mod- 
ulation and (b) for which x^ 7^ xi and a modulation 
appears, see figs. [2][b) and[3Ub). 

As the wavelength is further reduced, we find dramatic 
signatures of the presence of the FFLO phase, in both 
the cases (a) and (b) described above. Specifically, we 
find an unusually large response of the amplitude of the 
spin-dipole mode when the wavelength of the excitation 
becomes A = This response can be understood in 
terms of coupling to the gapless modes a.t k = = k* 
of the homogeneous system (see FiglT]). A small devia- 
tion from k* (of the order the inverse system size) allows 
mixing of this mode to the spin-dipole mode, causing the 
response at k* to be apparent in the dipolar motion of 
the atomic cloud. 

While the calculations we have presented are within 
RPA mean-field theory, as we now argue, the qualitative 
conclusions are valid more generally. Essentially, our re- 
sults rely on the fact that in the unpolarised FFLO phase, 
there exists low-frequency (of order luq) response that is 
sharply peaked at a spatial wavevector k*. In the above, 
this response was accounted for in terms of the broken 
translational symmetry of the mean-field state. As is 
well known, in a true ID quantum system no contin- 
uous symmetries are broken. Thus, the broken (phase 
and translational) symmetries of the mean-field FFLO 
state can lead only to power-law decay of the respec- 
tive correlation functions[l9|. However, the qualitative 
features described above are the same as those expected 
from an exact treatment of the system. The transition 
from the (unpolarised) superconducting phase to the par- 
tially polarised phase is marked by the closing of the 
spin-gap, leading to a second gapless sound mode. This 
can be viewed as a Luttingcr liquid representing the ex- 
cess fermionsfioj. These excess particles have density 
Pi — Pi, and thus a Fermi wavevector {kp^ — kp^). Thus, 
there are two gapless collective modes, arising from the 
fully paired particles, and the liquid of excess majority 
spin particles. Furthermore, as in the general theory of 
Luttinger liquids [2l| the spectral function of these ex- 
cess majority spin particles will show gapless response 
at multiples of twice their Fermi wavevector. This is 
2{kp^ — kp^) = 47r/AA which is precisely the wavevector 
A:* at which one finds the gapless response in mean field 
theory. Thus, the qualitative features of the collective ex- 
citation spectrum obtained in mean-field theory are fully 
consistent with those expected for the exact system. 

Any distinction between the exact results and those of 
mean-field theory will be further reduced in the finite-size 
systems on which experiments can be performed, where 
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the distinction between long-range and power-law cor- 
relations becomes blurred. Indeed, one can expect that 
density inhomogeneity will play a more significant role 
than will the power-law decay of correlations. Since the 
densities of spin-up and spin-down particles are inhomo- 
geneous, so too is (the local value of) Aa, so one can 
expect smearing of the condition A = Aa/2. However, 
even for the relatively small systems studied in Figs. [2] 
and [SI the collective mode response has a sharp onset 
at the value A = Aa/2 at the centre of the trap. The 
appearance, in Fig. [3I^b), of a large response at wave- 
lengths slightly smaller than Aa/2 is associated with the 
non-constancy of Aa, which decreases slightly towards 
the edge of the FFLO phase. 

We propose that the signatures of Figs. [2] and [3] pro- 
vide a convenient way to detect the FFLO phase in ex- 
periment. The sharp feature at Aa/2 appears only when 
the gap displays oscillatory behaviour, so is absent for a 
non-interacting gas and for temperatures when the FFLO 
phase disappears (fc^T ~ O.lEp [H]). 

In order to observe these signatures in experiments one 
needs to be able to create a variable wavelength optical 
lattice. This can be done using the technique used by 
Steinhauer et al. [i^l . While the signal appears for a spin- 
symmetric perturbation, the response is larger for spin- 
asymmetric. Thus, any spin-dependence of the optical 
lattice will enhance the ability to distinguish the peaks 
associated with the oscillation of A from the background 
peaks. Optical lattices with spin-dependence have been 
created by Mandel et al. [2J|, and would be ideally suited 



to this purpose. There are two natural ways to find the 
response of the spin-dipolc mode. One way is to fol- 
low precisely the approach of the calculation, and make 
the strength of the optical lattice time dependent, as in 
Ref . [is! . One can then selectively excite the spin-dipolc 
mode by bringing the temporal oscillation of the lattice in 
resonance with the spin-dipolc mode. Another way (sim- 
ilar to the approach used in Ref. [2^) is to apply a static 
periodic potential to the system, allow it to equilibrate, 
and then switch off this potential abruptly. This will 
excite collective modes of many frequencies from which 
the response of the spin-dipolc mode can be obtained by 
Fourier transform. 

In either case, we emphasise that the required mea- 
surements of particle density are on the length scale of 
the size of the cloud . Thus, our proposed method allows 
the probing of the microscopic physics by perturbing the 
system with short wavelengths while requiring only the 
measurement of densities on the length scale of the cloud. 
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